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. Using the methods introduced in a previous paper, we consider the dissociation of J/ip by the 

' lowest-lying pseudoscalar and vector mesons in the hadronic fireball formed in heavy ion collisions, 

^SJ ' assumed to be a hadron gas at temperature T. Absorption by nuclear matter is accounted for as well. 

We compare with the S-U and Pb-Pb data presented by the NA50 Collaboration. From data at low 
centrality we find T = 165 — 185 MeV, close to the predicted temperature of the transition to quark- 
' gluon plasma and to the temperatures measured by hadron abundances. We extrapolate to higher 

centralities with the approximation of scaling the energy density of the fireball with the average 
04 ■ baryon density per unit transverse area. Using the energy-temperature relation of the hadron gas 

made by the same pseudoscalar and vector mesons, the fall off of J/ip production shown by the 
(N ■ NA50 data can be marginally reproduced only for the highest temperature, T = 185 MeV. If we 

use the energy-temperature relation of a Hagedorn gas with limiting temperature Th = 177 MeV, 
predictions fall short from reproducing the data. These results suggest that a different mechanism 
is responsible for the J /iji suppression at large centralities, which could very well be the formation 
of quark-gluon plasma. 
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I , I. INTRODUCTION 

CD , „ 

In a recent paper [l| we have presented a theoretical study of the dissociation cross-section tt -t- J/t/" — > D^*'^ + £>(*^ 
to analyze the absorption of J/tj} in heavy ion collisions. The disappearance of the J ftp's created in the early stages 
of the collision was considered as mainly due to two subsequent processes: 
• 1. The interactions of J/ip with the nuclear medium traversed during the interpenetration of the two heavy nu- 
^ , clei,with an absorption cross section per nucleon experimentally determined 2] from the inclusive J/ip production in 
- - « p + A collisions; 

2. The dissociation process: tt + J/tp _D'*) -|- Z)'*^ induced by the pions of the hadron gas formed in the heavy-ion 
collision, the so-called co-moving particles 3]. 

In , j the absorption curves thus obtained have been compared to a data analysis by the NA50 collaboration 0]. 
Fitting the data at low centrality we found a value for the inner temperature of the pion gas of about T = 225 ± 
15 MeV. Extrapolation of the absorption curve at large centrality was done by assuming an increase of the energy 
density deposited in the fireball proportional to the average number of nucleons per unit of overlapping area Q , a 
quantity which increases with decreasing impact parameter. The results did not seem to reproduce the fast increase of 
absorption observed by NA50, thus lending some support to the idea that in large centrality collisions a quark-gluon 
plasma phase is produced, which would impair the formation of the J/ip because of QCD Debye screening j5|.As 
already noted in |l|, however, for this method of analysis to be accurate one should go beyond the approximation of 
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a simple pion gas and consider the more realistic case in which further particles/resonances appear in the fireball at 
thermal equilibrium, which can appreciably contribute to J/ij} dissociation. 

In the present paper, we extend our previous work in several directions. First, we include the J /ip dissociation cross 
sections by the lowest lying pseudoscalar and vector mesons (Sect. II). We use the Constituent Quark Model (CQM) Q 
as before, with couplings computed with flavour SU(3) symmetry, and nonet symmetry for the vector mesons. Real 
particle masses are used for reaction thresholds and to compute particle abundances in the hadron gas. 

After defining the different lengths that characterize the collision (Sect. Ill), the absorption lengths vs. temperature 
are computed in the heath bath made by the hadron gas of pseudoscalar and vector mesons at a given temperature 
T (Sect. IV). Not unexpectedly 0, the largest correction to the pion gas situation is due to the vector mesons. 
Notwithstanding the considerably larger mass, they contribute to J /if) absorption about as much as pions, for two 
reasons: (i) dissociation reactions have very low or no threshold, thus making all particles useful as opposed to pions, 
which are effective only above a relatively high-energy threshold; (ii) the large spin and flavour multiplicities (a factor 
of 24 with respect to 3 for pions, for the complex p-\- uj + K* + K*). 

The total inverse absorption length, 'Si{piai)T, is considerably increased with respect to the pure pion gas. Cor- 
respondingly, we find a substantial decrease in the estimated temperature of the fireball in low centrality collisions, 
which turns out to be now in the range 165 — 185 MeV (see Fig.[7J. Extrapolation to higher centrality is done using 
the relation between energy density and temperature appropriate to the hadron gas we are considering and it leads 
to increasing temperatures with increasing centrality up to 180 — 200 MeV at zero impact parameter, I = 2R — 13 fm 
for Pb {I is the linear size of the fireball and R the nuclear radius). The corresponding hardening of the attenuation 
curves, however, do not really reproduce the falling of J/ip production in this region, as shown by Fig.[7| supporting 
the idea of a change of regime, this time in a more realistic region of the thermodinamical parameters with respect 
to Q. 

In the range of temperatures we have found, one expects to deal with a hadron gas of increasing complexity 
approaching the Hagedorn gas , with an infinite number of resonances and a level density exponentially increasing 
with mass. This situation is known to give rise to a limiting temperature j^, which was interpreted in as the 
temperature at which the transition to a quark-gluon plasma phase starts to take place. In Sect. V, starting from 
the fit at low centrality, we extrapolate to higher centrality with the energy-temperature relation appropriate to the 
Hagedorn gas, with a Hagedorn temperature Th = 177 MeV. This temperature is consistent with the transition 
temperature estimated in QCD lattice calculations ^Jli with the empirical density of hadron levels up to 2 GeV, 
as estimated by 0, and with the temperatures determined at SPS and RHIC from particles relative abundances at 
freeze-out . We find that the increase of the energy deposited in the collision does not lead anymore to an increase 
of temperature, hence of the opacity: the extra energy deposited in the fireball for increasing centrality goes into an 
increase of the thermodinamical degrees of freedom, which ultimately should lead to the phase transition. 

The curve corresponding to the absorption by the Hagedorn gas falls definitely short to explain the further decrease 
of J/ijj abundance, as observed by NA50. It is tempting to interpret the discrepancy shown in Fig.|Slas a signal of 
quark-gluon plasma formation. However, caution is required, due to the unavoidable truncation we have made in the 
calculation of the opacities. While higher resonances are expected to contribute less and less, their cumulative effect 
could resum to a large contribution to the dissociation cross-section and lead to an increase in opacity even for the 
very small temperature increases allowed by the vicinity of the limiting temperature. 

Sect. VI contains conclusions and outlook. We find it very satisfactory that a microscopic calculation based on the 
CQM and the observed opacities produce values of the temperature that are consistent with the measured freeze- 
out temperature in ion collisions and are just at the border of the transition temperature predicted by QCD lattice 
calculation and by the Hagedorn hadron gas. A point borne out clearly by our analysis is that the picture would 
be considerably more clear if the relative normalization of S-U vs. Pb-Pb data was better understood this 
respect, further measurements at low energy, aimed at resolving the experimental issue of the relative normalization, 
could be very useful. 

All in all, the results shown in Figs. [3 (a), (b) or, even more, in Fig. [S] are quite impressive. They seem to us very 
suggestive that a new mechanism to suppress J/ip is setting in at large centralities. This could very well be the 
formation of quark gluon plasma. It would be extremely important to correlate other signals to the present one, in 
order to get to a definite conclusion about the presence of a phase transition. 

In this paper we correct few mistakes which affected our analysis in 1]. The pion cross sections presented here 
supersede those in where we found a trivial numerical mistake; the average lenght traversed by the J/ip in the 
fireball was incorrectlyestimated to be (6/10)/ rather than (3/8)L All together these corrections would increase the 
temperature given in 0] by 15 — 20 MeV, leaving all conclusions unchanged. Finally, a misprint led to the factor 27r/3 
in Eq. (13) of Ref. [3, to be replaced by tt/2. 
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II. CROSS SECTION COMPUTATIONS 

The evaluation of the cross sections for the processes (tt, 77, /C, p, lj, if *, 0) + J/%jj proceeds through 

the computation of diagrams very similar to those listed in Such tree level diagrams involve effective tri- linear 
53 = (tt, r], K, p, or 33 = J/ and four-hnear 54 = (tt, ry, K, p, ...) J/V'-D(*j-D(*j couphngs which we 

can estimate using the CQM model, originally devised to compute exclusive heavy-light meson decays and tested on 
a quite large number of such processes 

CQM is based on an effective Lagrangian which incorporates the heavy quark spin-flavor symmetries and the chiral 
symmetry in the light sector. In particular, it contains effective vertices between a heavy meson and its constituent 
quarks (see the vertices in l.h.s. of Fig. whose emergence has been shown to occur when applying bosonization 
techniques to Nambu-Jona-Lasinio interaction terms of heavy and light quark fields p^ . On this basis we believe 
that CQM is a more solid approach to the computation of 53, 174 if compared to the various methods available in the 
literature, often based on SU(4) symmetry (for a review see e.g. [TH). 
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FIG. 1: Basic diagrammatic equation to compute 33 and couplings. 



In Fig. ^ we show the typical equation which has to be solved in order to obtain .94(53) in the various cases at hand: 
on the r.h.s. we represent the effective four-linear couplings to be used in the cross section calculation (to obtain the 
tri-linear coupling we suppress either the J/ip line or the dashed line of the light particles); the effective interaction 
at the meson level (r.h.s.) is modeled as an interaction at the quark-meson level (l.h.s. of Fig. 

The J/V' is introduced using a Vector Meson Dominance (VMD) Ansatz: in the effective loop on the l.h.s. of Fig.Q] 
we have a vector current insertion on the heavy quark line c while on the r.h.s. the JJjb is assumed to dominate the 
tower of J^^ = 1 , cc states mixing with the vector current (for more details see [la|). Similarly, vector particles 
coupled to the light quark component of the heavy mesons p, w, when q — {u,d), or K* when one or both light 
quarks are q — s, are also taken into account using VMD arguments. The pion and other pseudoscalar fields have a 
derivative coupling to the light quarks of the Georgi-Manohar kind Tf] . 

Once established the form of the effective vertices occurring in the loop diagram on the l.h.s. of Fig. ^ one has 
just to compute it. The momenta running in the loop are limited by two Schwinger cut-off's: one in the ultraviolet 
and one in the infrared. The ultraviolet cut-off A has been set to the chiral expansion scale A^ ~ 47r/^ while the 
infrared fi prevents loop momenta to access the confinement energy region (CQM does not include any confining 
potential). To give a flavor of the kind of calculations involved, consider that the loop integral for the four- linear 
coupling pJ/t/jD'-*^ D^*") is written as follows: 



/ , rrz ^ f dH 

ZHrnHZn'mH' x / ,„ ,.. x 
J (2n)H 



Tr 



(1) 



where H and H' represent the heavy- light external meson fields labeled by their four- velocities v,v'] the Feynman rules 
for their couplings to constituent quarks, in this case the \/ZH'niHZH''mH' coupling factor, are discussed in The 
heavy quark propagators (those obtained by the standard Dirac propagator in the limit of very heavy Q = h,c mass) 
are also labeled by v, v' while we have the usual expression for the light propagators. The parameter A appearing in 
the heavy propagator is defined to be A = — ttiq, the mass of the heavy- light meson minus the mass of the heavy 
quark contained in it. A is the the main free parameter of the model. It varies in the range A = 0.3 — 0.5 GeV for 
u, d light quarks and 0.5 — 0.7 GeV for strange quarks [1^]. Varying A allows to estimate the theoretical error. 
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The p is coupled via its VMD coupling to the light quarks, e being its polarization. The J/i/", with polarization 77, 
is also coupled via VMD to the heavy quarks {rj appears in the trace between two heavy quark propagators, e appears 
between two light quark propagators). In front of this expression we have the fermion loop factor. 

The way to regularize and compute such integrals is discussed in Il9| . The trace computation in will introduce a 
number of scalar combinations of the momenta and polarizations of the external particles. Each of these combinations 
will be weighted by a scalar integral which amounts to a numerical factor: what we call the coupling. Actually such 
scalar integrals depend on the energy of the p. In general the expressions obtained in the <Jpj/^ computation appear 
to be quite complicated functions of Ep if compared to those obtained when studying only J interactions with 
pions |lg |. On the other hand, having in mind a hadron gas at a temperature T k, 170 MeV, we restrict our study to 
the low energy interval Ep ~ 770 — 1000 MeV and rely on the fact that the Boltzmann factor exp(— i5/T) will cut off 
the high energy tails and work as a high energy form factor in the thermal averages {p<j)t (we have explicitly tested 
that this occurs). 

The couplings pD^*^ D^*\ are computed in the same framework by simply dropping the J/tp line in the diagrammatic 
equation of Fig. ^ as mentioned above. 

The Lu particle contribution is introduced in the thermal average by simply multiplying the p cross-section by a 
factor of 4/3 (nonet symmetry). The amplitudes for K* are deduced from those of the p by SU(3) symmetry. However, 
computing the K* cross-sections, where a meson will be produced in the final state, we take into account exact 
masses and thresholds. Similarly rj and K contributions are computed by SU(3) symmetry from the pion couplings 
but accounting for different kinematical thresholds. The calculation of the (f) contribution is performed by replacing 
the constituent light quark mass, m — 300 MeV for q = {u, d), with a constituent strange quark mass m = 500 MeV. 
Such an operation requires accordingly the modification of the Zh couplings (see Eq. |^) and of the infrared cutoff, 
as described in [l8l |. 

The final results for cTttj/^ and Upjj,^ are displayed in Fig. [3 We give the cross sections for 7r+ and p+. 
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FIG. 2: The cross sections for the processes (7r+,p+) -|- J/-0 — > D(*'-D(*' versus energy. 



III. THE GEOMETRY OF THE COLLISIONS 

The geometry of the heavy ion collision is shown schematically in Fig. which depicts the time-evolution in the 
center of mass frame. 

The impact parameter, 6, is defined, as usual, as the transverse distance of the centers of the two nuclei. We 
consider the J/i/" to be created with Feynman's a; « 0, during the overlap of the two nuclei. These particles have to 
overcome absorption from the column density of nucleons. In the center of mass frame the length of the column is 
L/7. In the same frame, the density of nucleons is Pnuci. • 7, so that the absorption factor is Lorentz invariant and 
given by exp(— puuci.CTnuci.i)- This is the same formula used in Ref. Cnuci.j the nuclear absorption cross-section, 
has been determined by NA50 from the behaviour of the cross-section of p + A ^ J/^ + anything as function of A. 
We take from Ref. 0: 

CTnuci. = 4.3 ± 0.3 mb (2) 
Pnuci. = 0.17 fm-3. (3) 



(a) 




(b) 

A 



(c) 



L/y 




time 

FIG. 3: Time evolution of the collision of two heavy nuclei seen in the c.o.m. frame. A fireball with baryon number Nb = is left behind 
the two receding nuclei. The fireball is expected to thermalize at the temperature T. The nuclcons belonging to the nuclear overlapping 
regions fly apart as unbounded particles. The dimensions L/7 and I are respectively the linear dimension of the nuclear column the J/i/> 
has to traverse during the interpenetration of the two nuclei and the linear size of the fireball. 



After collision, nuclear matter appears in part as a cloud of free (wounded) nucleons from the initially overlapping 
parts of the nuclei, represented as dotted regions in Fig. O (c), in part as forward going fragments from the non 
overlapping regions. The NA50 Collaboration has installed a Zero-Degree-Calorimeter (ZDC) which measures the 
energy of the forward going fragments, proportional to their nuclcon number and therefore to their size, thereby 
determining the value of the impact parameter, 6, for each collision 4]. The relation between L and b has been given 
in 01 in the framework of the Glauber theory. 

In Fig.|31we show the hadron fireball produced by the central collisions of the interacting nucleons ^ (the comoving 
particles). The fireball has a transverse dimension, approximately equal to the length of the overlapping region: 

l = 2R-b, (4) 

where R is the nuclear radius. If we take, for simplicity, a spherical fireball, the average length that a J/tp has to 
traverse before leaving it is (3/8)/, so that the attenuation factor due to absorption by the comoving particles is: 



^comoving 



(X exp 



o 



(5) 



the subscript i labels the species of hadrons making up the fireball, pi the number density of the effective (i.e. above 
threshold) particles and ai the corresponding J/ip dissociation cross-section. Brackets indicate an average over the 
ener gy d istribution in the fireball. As noted before, we can express the nuclear absorption length, L, as a function 
of b [2y| and therefore, using I^J, as a function of I. Putting all together, we write the attenuation of the J/-0 as a 
function of I according to: 



N X exp[-pnucl.O-nucl.i(0] X ^Xp 



(6) 



where N is an appropriate normalization constant. As discussed in Rcf. 3], the energy density deposited in the 
fireball is proportional to the average flux of nucleons participating to the collision, i.e., the number of nucleons per 
unit transverse area. This quantity increases with decreasing impact parameter. A simple estimate of the ratio of the 
energy density for two different values of b was given in according to: 



E{b2/R) 9{b2/R) 



e{b,/R) g{b,/R)' 



(7) 



with the geometrical factor g{b/R) given by 



9{b/R) = 



(l-5/2i?)2(l + 6/4i?) 



2 arccos(V2i?) - {b / 2R) ^fT^WJlR^ 



(8) 
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IV. THERMAL AVERAGES AND ABSORPTION IN THE RESONANCE GAS 

The fireball depicted in Fig.|2|(c) is expected to quickly thcrmalize (see 0), giving rise, at least for low centrality 
collisions, to a hadron gas at temperature T. The thermalisation hypothesis is supported by the data from SPS 
and RHIC which show abundances and momenta distributions at freeze-out compatible with a temperature T = 
170 — 180 MeV, see e.g. Ref. and Ref. for a review of recent data. 

In Ref. we have assumed the thermalised fireball to be a pion gas. We consider now a more general case of 
a boson gas made by the lowest-lying pseudoscalar and vector mesons, retaining, for simplicity, the assumption of 
vanishing chemical potential. What matters, for our considerations, are particle densities above the dissociation 
threshold. Assuming the J/ij; to be at rest, this means particles with energy larger than Eth.' 



th. 



{Mf, 



2M 



(9) 



m is the projectile particle mass and Mf^^ the sum of the final particles masses. Explicitly, we have: 



PiT) 



N 
2^ 



dE- 



pE 



Eth.. 



■E/kT _ I ■ 



(10) 



N is the total multiplicity (spin times charge, TV = 3, 9 for pion and for p, respectively) and p = y/E^ — m? . We 
report in Fig. ^ and in Tables I and II the number densities above threshold for the particles considered. In Table I, 
we give in parenthesis the total number density of pions as well. 



Number density (fm"^) 
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0.08 
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0.04 
0.02 


ICQ 120 140 160 180 200 
T(MeV) 

FIG. 4: Number densities of the hadron gas particles h being above the kinematical threshold for the hj/ip — » d|*jD|*j reaction. 



The energy density associated to the various particles are listed in Tables III to V. In this case the sum over the 
full energy range is understood: 



The contribution of the other particles is not at all negligible with respect to pions, even at temperatures as low as 
T ~ 150 MeV, 9], particularly in terms of the energy density. 

We report in Figs. |S1 (a,b) the energy densities divided by eo = T'^tt^/SO. The ratio, g, gives the effective number 
of degrees of freedom which are active at the particular temperature. Pions give g = 3 already at T = 150 MeV 
(Fig. |S1 (a)), which is increased to g = 7 — 10 by the other particles, for T — 150 — 180 MeV. In Fig. [Sl(b) we specify 
the contributions of the different particles, with Kaons {N — 4) and the lowest lying vector mesons, p and oj {N = 12) 
providing the dominant contributions. 
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FIG. 5: The effective number of degrees of freedom g active at a particular temperature. In panel (a) we show the dramatic increase 
of g for the Hagedorn gas with limiting temperature Tn = 177 MeV, see 1151 . In panel (b) we specify the contribution of the different 
components of the hadron gas. 



Thermal averages of the product p ■ a give the inverse of the absorption length due to each particle species, x, and 
are done with analogous formulae: 



, N pEa{E) 

{p-'^x+.J/i,^DMDM)T - ^ dE-^jj^^^r—j. (12) 

We give the values of the inverse absorption lengths as function of temperatures in Table V and Fig. p and oj 
give larger contributions than the pions. This is due to the absence of threshold in the dissociation reaction, which 
makes particles of all energies to be effective, and to the large multiplicity. 




150 160 170 180 190 200 210 220 



T(MeV) 



FIG. 6: The inverse absorption lengths as a function of temperature. 



The absorption lengths thus computed are inserted in the absorption master formula @ . We compare in Fig. E the 
absorption curves thus obtained with the NA50 data in S-U and Pb-Pb coUisions, for T = 165 - 185 MeV |2l|. The 
normalization constant in © is chosen to fit the data around I = 3 fm. The data for ^ < 4 fm favour temperatures 
of this order. In the same figure, we show the curve corresponding to the pure nuclear absorption. 

To extrapolate to higher centralities, we keep into account the increase in the energy density deposited in the 
collision as explained in Sect. III. The increase in temperature for increasing centrality makes the absorption to 
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increase, as shown in Fig. [71(b), with respect to Fig. |7|(a) where the geometrical effect is neglected. In Fig. |7|(b) we 

label the curves with the value of the temperature at I — 3.4 fm; the temperature increases along the curves, up to 
T w 180 - 200 MeV for I = 11 fm. 




2 4 6 8 10 12 2 4 6 8 10 12 
l(fm) l(fm) 



FIG. 7: The absorption curves are compared to the most recent data analysis of the J/i/i yield in S-U (star) and Pb-Pb (box) collisions ll.*^ . 
Left panel: the absorption by comoving particles is described by a single temperature; the effect of nuclear absorption alone is shown 
by the upper curve (since the relation between L and I is non-linear, the nuclear absorption does not give rise to a simple exponential, 
unlike the isothermal absorption by comoving particles). Right panel: the temperature increase due to the geometrical effect described in 
Sect. Ill is accounted for with the energy-temperature relation of the hadron gas. The temperatures indicated refer to I = 3.4 fm. 

The inclusion of the lowest lying resonances improves greatly the consistency of the whole picture, bringing the 
estimated temperatures below the temperature expected for the phase transition and closer to the temperatures 
measured at freeze-out, from the abundances of the different hadrons. The curve labeled with T — 175 MeV, which 
has T varying from 175 to about 190 MeV, fits well the data for low centrality but still falls short from reproducing 
the observed drop in the production of the J/^ above, say, I = 5 fm. The curve with initial temperature T — 185 MeV 
fits the low centrality data as well and follows closer the data at large centrality. However the temperature exceeds 
200 MeV at / ~ 11 fm which is likely too a high value for a hadron gas (see next Section). 

The increase in temperature due to the increase in energy density that we find is definitely less pronounced than 
what found for the pion gas in Q . The reason is that in the resonance gas the number of degrees of freedom increases 
appreciably with temperature. The extra energy density provided has to be shared among more and more degrees of 
freedom and the temperature increases less than with a fixed e = CT* power law. This behaviour begins to reproduce 
what is expected in the case of a Hagedorn gas, with an exponentially increasing density of resonances per unit interval 
of mass which we turn now to consider. 



r (MeV) 


(p^rot) (fm-^) 




p« (fm-^) 


(fm-^) 


{pa}^ (fm-^) 


150 


0.004 (0.12) 


0.001 


0.006 


0.004 


0.001 


165 


0.008 (0.17) 


0.003 


0.012 


0.009 


0.003 


180 


0.016 (0.22) 


0.006 


0.022 


0.017 


0.007 


195 


0.028 (0.29) 


0.010 


0.039 


0.030 


0.012 


210 


0.045 (0.37) 


0.016 


0.062 


0.050 


0.021 


225 


0.069 (0.47) 


0.024 


0.095 


0.080 


0.034 


240 


0.102 (0.57) 


0.035 


0.139 


0.121 


0.052 



TABLE L Number densities and inverse absorption lengths for the pseudoscalar particles in the gas. The number densities 
are relative to those particles which are over the kinematical threshold to open the _D|*jI)|*j channel. In the case of pions, we 
report in parenthesis their total number density (i.e., the one including pions below threshold). 
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T (MeV) 








{paY+- (fm-i) 




(pa>^ (fm-i) 


150 


0.018 


0.017 


0.001 


0.007 


0.002 


0.0007 


165 


0.037 


0.035 


0.003 


0.016 


0.006 


0.002 


1 OA 

ioO 


U.0d8 


0.064 


0.00/ 


n AO K 

0.030 


A A1 A 
0.014 


A AAK 

O.OOo 


195 


0.115 


0.108 


0.012 


0.070 


0.030 


0.010 


210 


0.184 


0.172 


0.021 


0.131 


0.059 


0.020 


225 


0.279 


0.260 


0.034 


0.231 


0.108 


0.039 


240 


0.404 


0.377 


0.052 


0.388 


0.188 


0.070 



TABLE II; Same as Table I for vector particles. 



T (MeV) 


e'^ (MeV/fm-^) 


e" (MeV/fm-^) 


(MeV/fm^'O 


150 


58.8 


6.6 


31.4 


165 


87.5 


11.5 


53.5 


180 


125.5 


18.8 


85.3 


195 


175 


28.9 


129 


210 


237 


42.6 


188 


225 


314 


60.6 


265 


240 


409 


83.6 


362 



TABLE III: Energy density for the pseudoscalar particles in the gas. 



V. ABSORPTION IN A HAGEDORN GAS 



In the range of temperatures we have found, one expects to deal with a hadron gas of increasing complexity, 
approaching the Hagedorn gas 8], which has an infinite number of resonances with a level density exponentially 
increasing with mass. This situation gives rise to a limiting temperature 0, the Hagedorn temperature Tjj, which 
was interpreted in (1(TI | as the temperature at which the transition to a quark-gluon plasma phase starts to take place. 
The Hagedorn temperature can be estimated from a fit to the resonance level density below 2 GeV, appropriately 
smeared down in the lowest end. In |^ a value Th — 158 MeV is estimated. However, this value is lower than 
the observed hadron temperature at freeze-out and also lower than the estimates of the transition temperature from 
lattice QCD calculations, see Ref. [m. As a reasonable compromise, we take = 177 MeV, which still fits the 
resonance spectrum and it agrees with lattice QCD calculations and with the observed freeze-out temperature. We 
use the partition function of the Hagedorn gas in the form |^: 

with Hi: 



C = 2.12GeV^; mo = 0.96 GeV; T/^ = 177 MeV (14) 
and the energy density {(3 — 1 /T) : 

e = -^\nZH + eAT). (15) 

To get the total energy density we have added the contribution of the pion, Eq. (|ll|l which is not included in p{m). 
The number of effective degrees of freedom vs. temperature is shown in Fig. 

To make a numerical calculation possible, we have to make an assumption on the J /i}) dissociation cross-section by 
the hadron resonances. We assume that only the pseudoscalar and the vector mesons we have considered before are 
relevant to dissociate the J/V'- Starting from the results of Fig. [7| we extrapolate to increasing centrality using the 
energy-temperature relation of the Hagedorn gas. This is shown in Fig.|Sl using as initial temperature T = 175 MeV. 
The result is quite spectacular. The sharp rise of the degrees of freedom due to the vicinity of the Hagedorn temperature 
makes so that the temperature of the gas practically does not rise at all, the dissociation curve cannot become harder, 
and the prediction falls really short from explaining the drop observed by NA50. The simplest interpretation of 
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T (MeV) 


6"+'^ (MeV/fm^^) 


e^" {MeY/im-'-^) 


e'* (MeV/fm-^) 


150 


33.4 


31.4 


2.7 


165 


65.2 


53.5 


6.1 


1 on 


lib 


00.3 


loo 


195 


193 


129 


22.4 


210 


304 


188 


38.0 


225 


455 


265 


61.0 


240 


657 


362 


93.4 



TABLE IV; Energy density for the vector particles in the gas. 



T (MeV) 


p''°' fm-^ 


{poYr"' (fm-^) 


e"'°* (MeV/fm-^) 


g 


150 


0.050 


0.016 


153 


7.3 


165 


0.1 


0.037 


265 


8.7 


180 


0.18 


0.078 


436 


10.0 


195 


0.32 


0.153 


684 


11.4 


210 


0.51 


0.282 


1030 


12.9 


225 


0.77 


0.492 


1497 


14.2 


240 


1.13 


0.821 


2112 


15.5 



TABLE V: Total number densities, inverse absorption lengths, energy densities and number of active degrees of freedom, 
summed over the various components of the hadron gas. 



Fig. |S1 is that with increasing centrality, more energy is deposited but this goes into the excitation of more and more 
thermodinamical degrees of freedom leading to the final transition in the quark-gluon plasma. The curve shown in 
the figure would represent the limiting absorption from a hadron gas, anything harder being due to the dissociation 
of the J /if) in the quark-gluon plasma phase. 

Some words of caution are in order. In the framework of the CQM, it is certainly reasonable to expect the relevant 
insertions in the quark loop of Fig.^to correspond to the Dirac matrices 5", P, A, V, T and the latter to be dominated 
by the lowest qq, S'-wave states we have been considering. On the other hand, we cannot exclude that decreasing 
couplings of the higher resonances may eventually resum up to a significant effect, which would change the picture 
obtained by truncating the cross section to the lowest levels. 

However, in all cases where this happens, like e.g. in deep inelastic lepton-hadron scattering, the final result repro- 
duces what happens for free quarks and gluons. In our case, this would mean going over the Hagedorn temperature 
from the hadron into the quark and gluon gas, which is precisely what Fig. |S1 seems to tell us. 



VI. CONCLUSIONS 



The prediction that its formation is suppressed in quark-gluon plasma (5j makes the J /ip a very interesting probe 
for a possible phase transition in heavy ion collisions. When the prediction was made, it was believed that J/t/j would 
suffer negligible absorption in the nuclear matter and in the hadron fireball formed during the collision. Subsequent 
studies showed that this is not the case and, for the probe to be efficient, one must determine in a reliable way the 
absorption of J /ip from these two "conventional" sources. The absorption length of J /ip in nuclear matter has been 
determined |2| from the the cross-section oi p + A-^ J/ '(/'-(-anything as a function of the nucleon number, A. 

A method to estimate the absorption from the comoving particles was presented in a previous paper 0], under 
the simplified assumption that the comoving particles are a thermalised pion gas at temperature T . Starting from 
a calculation of the tt -|- J/i/i £)(*)/)(*) cross-section based on the Constituent Quark Model |^, we showed that 
the absorption is potentially large and strongly temperature dependent. T itself can be determined by a fit to the 
data at low centrality, where the pion gas approximation could be better justified. Comparing the extrapolation to 
larger centralities with data, one could then judge if a further suppression is at work. The temperatures found in Q], 
T Ri 200 MeV, were rather on the high side but still in the right ball-park, a reassuring indication for a calculation 
completely based on microscopic parameters. On the other hand, for the method to be reliable, it is necessary to 
apply a better approximation than the simple pion gas. 

In the present paper we have extended the calculation to the lowest-lying pseudoscalar and vector mesons, the 
S'-wave qq states. Contributions from higher resonances like the Ai and the other P— wave qq states could also be 
considered, at the price of increasing theoretical uncertainties, and should not change the picture. Vector mesons give 
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FIG. 8: The same as Fig. |7| without (a) or with (b) geometrical effects taken into account, using the energy-temperature dependence of 
the Hagedorn gas. The sharp rise of the degrees of freedom due to the vicinity of the Hagedorn temperature makes so that the temperature 
of the gas practically does not rise at all, the dissociation curve cannot become harder, and the prediction falls short from explaining the 
drop observed by NA50. Nuclear absorption alone is indicated by the upper curve. 



substantial contributions to J/i/j dissociation, due to the higher multipHcities and because they are close or above 
threshold. As a consequence, we found a larger absorption in the heath bath at temperature T and a much lower range 
of temperatures to fit low centrality data, T = 165 — 185 MeV. This range agrees with the temperatures found from 
(i) particle abundances in ion collisions at freeze-out, Qi) the transition temperature found in lattice QCD calculations 
and (iii) the limiting temperature of a hadron gas |a ISl i estimated from the experimental hadron level density. 

In the range of temperatures considered, the number of active degrees of freedom is rapidly increasing, due to the 
excitation of more and more hadron levels (Fig. [S]) . This concept provides the ground for the existence of a limiting 
temperature in a hadron gas, the Hagedorn temperature. For the gas made by the comoving particles, the increase of 
the degrees of freedom "blocks" the temperature to the values observed for the low centrality events and prevents the 
opacity of the fireball to J/t/i to increase further. In fact, assuming a limiting temperature of about 170 — 180 MeV 
from the outset, the decrease of J/'0 production observed by NA50 at large centralities cannot be explained with 
absorption by either the nuclear matter or the comoving particles, as shown in Fig. |S1 

We cannot exclude that an infinite tower of higher resonances with decreasing couplings could eventually change 
the picture obtained by truncating the cross section to the lowest levels. Should this happen, however, we would 
expect the result to be close to that obtained by replacing the resonances with quarks and gluons, i.e. the situation 
one encounters in the QGP. 

In conclusion, we believe that our calculation has produced a reliable estimate of J /ip absorption in a hadron gas 
where only the lowest resonances are excited. It can agree marginally with the data at large centrality, at the expenses 
of accepting temperatures which are quite high (T « 200 MeV) with respect to the current theoretical estimates of 
the limiting hadron temperatures. On the other hand, should one introduce the limiting temperature as a boundary 
condition, it is impossible to fit the data at large centrality, thus providing considerable support to the idea that the 
suppression is produced by the quark-gluon plasma. 

Further progress can be envisaged in several directions. Further measurements at low energy, aimed at resolving 
the experimental issue of the relative normalization of S-U vs. Pb-Pb data, would restrict the range of temperatures 
required at low centrality and allow for a more precise extrapolation with respect to what shown in Fig. [7| The 
onset of other possible signatures of the phase transition, such as strangeness production, should be correlated to the 
present signal, in centrality and energy. Finally, a quantitative calculation of the expected absorption curve, should 
the QGP be formed, would be of crucial importance. 
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